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Abstract 

A new linearization approach incorporated with chiral symmetry is 
established and applied to study the irK elastic scatterings. We demon- 
| strate that the k resonance exists, if the scattering length parameter in 

^ ' the 1=1/2, J=0 channel does not deviate much from its value predicted 

Q^ ■ by chiral perturbation theory. The mass and width of the k resonance 

is found to be M K = 594 ± 79MeV, T K = 724 ± 332MeV, obtained by 
fitting the LASS data up to 1430MeV. Better determination to the pole 
C^) • parameters is possible if the chiral predictions on scattering lengths are 

taken into account. 
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X. 

c3 ■ 1 Introduction 



There have been many works devoted to study the resonance structure in the 
/ = 1/2, J = channel of nK scatterings. It has been suggested a long 
time ago that there may exist a resonance named k in this channel pQ. Very 
recently, both the E791 Collaboration [2] and the BES Collaboration jH] have 
found evidences for the k resonance, in the D — > Kirn and J/^f — » K*Kn 
channels, respectively, which makes the topic even more interesting than ever. 
Previous theoretical studies are mainly based on the LASS data |3] and the 
SLAC data [Hj for ttA' scatterings, from which the k resonance was suggested 
to exist in the J=1/2,J = channel, based upon various model dependent 
analysis (see for example Refs. jO]- [0]), whereas contradicting opinions also 
exist (see for example Refs. [10] -[121), casting doubt on the existence of this 
resonance. 



1 Address after Sep. 1st, 2003: Department of Physics, McGill University, 3600 University 
Street, Montreal, Quebec. H3A 2T8, Canada. 
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The reason that different and sometimes contradicting conclusions exist is 
partly, if not mainly, due to the fact that, this channel, like the I = J = 
channel in tttt scatterings, is of strong interaction nature and undergos strong 
unitarity corrections. Indeed, ttK scattering amplitudes have been calculated 
up to 1-loop order in chiral perturbation theory (x?T) [T31 ll4|. However, since 
chiral loop expansion is an expansion in terms of external momenta and masses 
of pseudo-scalar mesons, the perturbation series is in principle useful only at 
very low energies due to the violation of unitarity in perturbation calculation. 

This paper is devoted to the study on the k resonance in the s wave irK 
scattering processes. 2 For this purpose, we first develop some new disper- 
sion techniques which improves and extends our previous results in this di- 
rection El EI]- The main improvement is that in our present scheme, 
unitarity is manifestly preserved. Different contributions from poles and cuts 
to the scattering phase shift are classified, and different contributions to the 
phase shift are additive. Then we make use of the xPT results to estimate the 
left hand cut in the un-physical region. In this region - since it is further away 
from those resonance poles - ;\TT is expected to work well, at least qualita- 
tively. We find that the k resonance exist, if the scattering length parameter 
in the 1=1/2 channel does not deviate much from its value predicted by chiral 
perturbation theory. 

This paper is organized as the following: sec. Q is the introduction. In 
sec. 121 we review some background knowledge being used in the later dis- 
cussions. They include an introduction to the kinematics for irK scatterings 
and ;\TT results at 1-loop level. Especially we discuss the dispersion tech- 
niques previously developed for studying 7T7t interactions [THl HZj and modify 
those dispersion relations to meet the new kinematics. In sec. |3] we make 
a pedagogical analysis to the Pade approximation to the xPT amplitudes. 
This method, and its variations |19| . have been widely used in the literature 
to study the non-perturbative dynamics involving pseudo-Goldstone bosons. 
However, we reveal severe problems this approximation method encounters in 
ttK scatterings, similar to what happens in tctt scatterings [201 ■ We conclude 
that in dealing with chiral amplitudes, in some cases the Pade approximation 
is a poor unitarization method. Following the idea proposed in Ref. ^H] in 
sec. E] we develop a new method of unitarization which respects all known fun- 
damental properties of S matrix theory, which are unitarity, analyticity and 
causality, though crossing symmetry is not implemented automatically. Also 
efforts have been made to combine chiral symmetry and the results from chiral 
perturbation theory in the new unitarization scheme. The new unitarization 
scheme starts from first principles and is formally rigorous. Of course, approx- 
imations have to be made once it is used in practice, but our scheme clearly 
shows where those systematic errors induced by approximations come from 
- a property many models do not maintain. In this section we also discuss 
one issue with respect to the Breit-Wigner description of resonances. Sec. E] 

2 Part of the results are already presented in Ref. [13] . 
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devotes to the numerical fit to the ttK scattering data. A major topic in this 
section is the estimation to the background contributions, which is essential, 
and even vital, in studying broad resonances. In the rest of sec. |S] we present 
detailed numerical analysis on the location of the k pole with or with out fur- 
ther constraints from chiral perturbation theory. Sec. El is for discussions and 
conclusions, including a comparison to several other related works found in the 
literature. 



2 Preliminaries 

In this section we briefly review some basic properties of nK scattering ampli- 
tudes. In Sec. 12.11 we introduce the kinematics for txK scatterings. In Sec. 12.21 
we briefly review the known results on t:K scatterings from 1-loop 577(3) chi- 
ral perturbation theory. In Sec. l2.3l we establish the dispersion relations for ttK 
scatterings following the method proposed in Ref. [THj, some of the contents 
are new. 

2.1 Kinematics for ttK Scatterings 

The center of mass momentum in the s-channel is written as, 




where 

s R = (m K + m 7T ) 2 ,s L = (m K -m^) 2 , (2) 
and p(s) is the kinematic factor: 

p(s) = = -y/(s-s R )(s-s L ) . (3) 

V s s 

The partial wave S matrix is, 

S(s) = l + 2ip(s)T(s) , (4) 

where and hereafter we omit the indices of isospin and spin, I and J, of the par- 
tial wave amplitudes unless it causes confusion. In the single channel physical 
region the S matrix is unitary which leads to, 

ImT( S ) = p(s)\T(s)\ 2 , {&S = 1). ( 5 ) 

Since S is unitary in the single channel physical region, the partial wave S 
matrix can be conveniently parameterized as the following, 

S(s) = e 2tSl{s) , (6) 
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Figure 1: The left hand cut, circular cut and the right hand cut of irK 
scatterings. 

where 8% is the ttK scattering phase shift. The subscript 1 indicates that 5i is 
only defined, at this moment, in the single channel physical region, i.e., above 
the nK threshold but below the Krj threshold where 5i is real. Since in the 
I = | channel the partial wave 5* matrix above the Kr\ threshold is no longer 
unitary, the S matrix is then parameterized as the following familiar form, 

S(s) = V (s)e 2iS ^ , (7) 

where 62 is the phase shift experimentally measured above the Krj threshold, 
and rj(s) is the inelasticity parameter. Apparently Si and S2 are different ana- 
lytic functions. In fact they have a simple relation above the second threshold, 

S^s) = S 2 (s) - ~ log7?(s) , or 5 2 (s) = ReSi(s) (s > (M K + M v ) 2 ) . (8) 

The relation Eq. (JHJ) will be useful in later discussions. For convenience we 
in the following will also often use a simpler notation 5 to describe the ex- 
perimentally observed phase shift regardless which region it is defined, as is 
usually adopted in the literature. 

Because 7rK scattering is an unequal mass scattering process, the singular- 
ity structure of its partial wave amplitude is more complicated than the equal 
mass scattering. The cut structure for the partial wave amplitudes is depicted 
in fig. □ |2U 

2.2 ttK scattering amplitudes in SU(3) chiral perturba- 
tion theory 

Since the 7r and K has isospin 1 and ^, respectively, there are two independent 
ttK scattering amplitudes: T* K and T* K . To be specific, consider the process 
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which has purely isospin |, 

7T + ( Pl ) + K + ( P2 ) -> 7T + (J9 3 ) + K + ( P4 ) . (9) 

3 

The amplitude T* K depends on the conventional Mandelstam variables 

s = (pi + p 2 ) 2 , * = (pi - P3) 2 , w = (pi - P4) 2 , (10) 

with the constraint s + t + u = 2(ml + m^). The s « — > u crossing of the 
amplitude Eq. © will generate the process n + K~ — > 7r + /\~~, which has both 
7 = | and I — | components, so we have the following relation: 

= |drM,*) - ^(M,") . (11) 

3 

The amplitude T* K in £77(3) x ££7(3) chiral perturbation theory was given in 

Refs. nans, 

T w V(s, t, u) = T 2 (s, t, u) + Tf(s, t, u) + T[(s, t, u) + Tf {s, t, u) , (12) 

where T 2 denotes the tree- level part, Tj denotes the tadpole terms of order 
77 4 , T[ denotes the polynomial terms of order E A and T% denotes the unitary 
corrections: 

T 2 (s,t,u) = 1 (ml + m 2 K - s), 

I 

7f(s,t,u) = (m 2 K -ml)(-3fi n + 2[i K + /i v ) 

T[(s,t,u) = J ^(AL 1 r (t-2ml)(t-2m 2 K ) 

+2L 2 r ((s -ml- m 2 K ) 2 + (u - ml - m 2 K ) 2 ) 
+L 3 r ((u -ml- m 2 K ) 2 + (t- 2ml)(t - 2m 2 K )) 
+AL A r {t{ml + m\) - 4m>^) 
+2L 5 r ml (ml - m\ - s) + 8(277; + L 8 r )mlm 2 K ) 

T?(s,t,u) = ^^(t(u-s)(2MUt) + M r KK (t)) 

- t){L, K {u) + L Kv (u) - u{Ml K {u) + M r Kv (u))) 
+ (m 2 K -ml) 2 (M: K (u) + M r Kri (u))) 
+-(m 2 K - ml){K« K {u){hu - 2m\ - 2m 2 K ) 
+K Kri (u)(3u - 2m T 2 - 2m 2 K )) + J^ K (s)(s - m 2 n - m 2 K 

+i j; k (u)(uu 2 - i2«K + O + 4 « + O 2 ) 
+f - 1 W + O) 2 + ^(*)*(2« - ™ 2 ) 
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where the explicit expressions of functions fip and M P q, Lpq, Kpq, J P q are 
displayed in Ref. [221 • In the above Eq. ([T3]) of isospin I = | scattering 
amplitude, we rewrite the expressions of T 2 (s,t,u) and Tl(s,t,u) in terms of 
F 1 Fk rather than , following the conventional wisdom. Besides, two typos of 

the expression of T±(s, t, u) in Ref. [2] are corrected as being done in Ref. p. 
The partial wave expansion of the isospin amplitudes is written as 

T\ Sl t) = 16tt^(2J+ 1)Pj(cos 0)Tj(s) . (14) 
j 

The partial wave amplitudes in x?T expanded to 0(p 4 ) are 

Tj(s)=Tj t2 (s)+Tj A (s) . (15) 
The expressions of the partial wave amplitudes are very voluminous, and we 

2 I 

only give the tree level results of T 2 2 (s) and T 2 2 (s), 
T i fe \ _ rn K 2 + 711^-3 

I 0,2\ S ) — 



32F K F n 7T 

m h , x -3(m x 2 -m,, 2 ) 2 -2(m A - 2 + m 7r 2 )s + 5s 2 

T °' 2(S) = IT&FkF^s • (16) 

The above equations indicate a partial wave T matrix zero at s = m 2 K + m 2 
in the 1=3/2 channel and s ~ m 2 K — in the 1=1/2 channel. 3 At 1-loop 

level the locations of these Adler zeros will receive corrections depending on 
the Li parameters. These corrections are however rather small since the chiral 
expansion works rather well in the energy region. 



2.3 The Dispersion Representations for ttK Scattering 
Amplitudes 

Following the method of Refs. we define two functions F and F, 

F < s > - 2&r){w-im) (17) 

For isospin |, F and F have no right hand cut in the energy region we are 
concerning. For isospin |, F and F have the cut starting from (Mk + M^) 2 , 
but the right hand cut is very weak until the (Mk + M^) 2 threshold is reached. 

3 The second equation of Eq. 116(1 affords another tree level zero on the negative real axis. 
However after performing the partial wave integration, the amplitude receives an imaginary 
contribution at 1-loop level due to the left hand cut and the T matrix zero disappears. 
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The functions F, F define the analytic continuation of the phase shift 5 in the 
following way: 

cos(2«Ji(s)) = F(s) , 

sin^s)) = p(s)F(s) . (18) 

To evaluate various contributions to phase shifts from different singularities, we 
can construct dispersion relations for F(s) and F(s) like what have been done 
in Ref . fHl E! • However the dispersion relations for tcK amplitudes are more 
complicated than those for nir scatterings, since in here we have the circular 
cut at \s\ = M\ — M 2 as depicted in fig. [T] Considering only F for simplicity, 
we first take the contour of the left hand integral outside the circular cut, 

m = F(so)+ E -(.-yMl^ri-*-"!) W) 

s^tat [S ~ Sp){S - S p ) 7T J -co {Z-S)(Z-S ) 

(s — s ) F((m 2 K - ml + e)e id )(m 2 K - ml)e i9 d8 



2n Jn ((m 2 K — m\)e %e — s)((m 2 K — m 2 )e ie — s 
s-s ) f°° lm R F(z) 



7T J(m v +m K )2 (z — s)(z — S ) 



dz . (19) 



In the above relation, the integration on the circular cut is along the outer edge 
of the circle \s\ = m 2 K — ml and the sum runs over all poles, denoted as s p , 
outside the circle. If we consider the whole complex s plane, the contribution 
from the inner circular cut will be counteracted by the poles inside the circular 
cut. So we have the following relation: 

Q ~(s - s )Res[F(s p )} | (s - s ) /-fr*-™*) 2 Im L F(z) ^ 

sf^in (S - S P )(S - S p ) 7T J -(m^-ml ) (z - s) (z - 8 ) 

(s - go) r F((m 2 K -ml- e)e ie )(m 2 K - ml)e ie d6 

2tt J-n {{m 2 K - ml)e ie - s){{m 2 K - m\)e e - s ) ' 1 J 

where the sum runs over all poles inside the circle. Combining Eqs. (J19)) and 
(HD we get, 

s v all \ S - S p)\ S - S p ) 7T J-oo [Z - S){Z - S ) 



s ) 71 discF(z)dz (s - s ) r°° lm R F(z) ^ 



2iri (z — s)(z — so) n J(m v +m K ) 2 (z — s)(z — s ) 



Corresponding relations for F can be similarly written down. 

It is already possible to use Eq. (J2T|) to study the k pole problem, parallel 
to what is done in Ref. jTH] to study the a resonance. However in the following 
we will adopt another improved method which automatically incorporate the 
unitarity constraint. 



(21) 
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9 41 9fiS -1- 1 477 i 


-1 8 4Q98 -1- 1 661 1 i 


Jo 


Resonance 


0.759984 + 0.297399 i 


-0.144373 - 0.613777 i 


Resonance 


0.083424 + 0.449582 i (inside) 


0.00286773 - 0.0699666 i 


SPSR 


0.0610924 + 0.0188401 i (inside) 


0.0891571 + 0.673712 i 


SPSR 


0.299636 + 1.24856 i 


-2.18809 - 4.04587 i 


SPSR 


0.152309 + 0.343659 i (inside) 


0.0634636 - 0.133458 i 



Table 1: Resonances and spurious physical sheet resonances (SPSR) predicted 
by the [1,1] Pade approximant of ttK scattering on the complex s plane in both 
/ = | and I = \ channels using the values from Eq. ([24)1 . The pole position 
^= M + iT/2. 



3 The Unitarization of the Scattering Ampli- 
tudes — Pade Approximation 



Since the scattering amplitudes from chiral perturbation theory only satisfy 
unitarity perturbatively and can not predict the physical resonances by itself. 
We may restore unitarity by constructing the [1,1] Pade approximant for the 
partial wave amplitudes of ttK scattering, 

™M = TTISiT ' < 22 > 

1 T 2 (s) 

The [1,1] Pade approximant given above satisfies elastic unitarity if perturba- 
tive amplitudes satisfy the elastic unitarity order by order, 

ImT [1 ' 1] (s) = p(s)|T [1 ' 1] (s)| 2 . (23) 

We take one set of values of the low-energy constants L*(xl0 3 ) from Ref. PEj: 

L\ = 0.65 ± 0.28 , L 2 r = 1.89 ± 0.26 , L 3 r = -3.06 ± 0.92 , 
Lt r = 0.0 ±0.5 , L\ = 2.2 ±0.5 , 

L G r = 0.0 ± 0.3 , L s r = 1.1 ± 0.3 . (24) 

As an educative example, we use these values of low-energy constants to study 
the singularity structure of the S matrix of Pade approximation in the complex 
s plane. Poles of the S matrix in the first and second sheet in both I — | and 
I = \ channels are listed in table ^ Since the S matrix contains a circular 
cut at \s\ = {m 2 K — m 2 ), in the table we indicate explicitly whether the pole 
locates inside or outside the circular cut. From tabled we can find that the 
[1,1] Pade approximation predicts the existence of k resonance in / = ~ chan- 
nel, M + iT/2 = 0.759984 ±0.297399iGeV, but it also predicts many spurious 
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0.2 H . 1 1 1 1 1 1 1 1 1 1 1 H 1 . 1 . 1 ■ 1 1 1 1 1 1 1 

0.6 0.8 1.0 1.2 1.4 1.6 1.8 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

s 1/2 (GeV) s 1,2 (GeV) 

Figure 2: Various contributions to cos2<5q^ 2 and sin25o from spurious poles 
and left hand cut in the isospin | s wave Pade amplitudes. The Lj parameters 
are taken from the central value given in Eq. (J24j) . Here 'background' means 
all contributions on the r.h.s. of Eq. (}2Tj) except the pole terms. 



poles. The effects of these spurious poles are not always small. Especially 
in the I = | channel, the phase shift should have been given entirely by the 
background (subtraction constant + cut integrals) contributions, but in the 
Pade approximant it is dominantly contributed by the spurious pole contribu- 
tions (see the figure for cos25 3//2 in fig. |SJ). Various contributions of physical 
poles, spurious poles, left hand cut (including the circular cut) and right hand 
cut are clearly shown in figs. |2l and El Beside the problems mentioned above, 
Pade approximation also fails to give the correct s dependence at s — 0. 
Even though the related studies have made remarkable successes in, for ex- 
ample, correctly predicting various physical poles, [TH] we conclude that such a 
unitarization approximation contains apparent shortcomings, especially in the 
isospin 3/2 channel. 



4 A New Approach of Unitarization 

It is therefore necessary to find an alternative approach bridging correctly the 

5 matrix theory and perturbation theory. It will be the main purpose of this 
section. In here we generalize the discussion made in Ref. ^H] to the general 
case of unequal mass scattering and make a more complete analysis to the new 
unitarization scheme. 

4.1 Simple S Matrices 

Single channel unitarity of the S matrix tells us the following identity: 

cos 2 26(s) + sin 2 25(s) = 1 , (25) 
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Pade • 



0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 

s 1/2 (MeV) s 1/2 (MeV) 

1/2 1/2 

Figure 3: Various contributions to cos25 and sin 2% from the k pole, 
other spurious poles, left hand cut and right hand cut in isospin | S wave nK 
scattering. Here 'spurious poles' means all poles other than the k pole. 



which is the generalized unitarity relation and holds on the entire complex 
s plane. A physical S matrix is very complicated since it has various poles 
and cuts. We overleap this complicated S matrix and at this moment only 
consider some simple circumstances. The word "simple" means that the S 
matrix contains either only one pole or zero on the real axis, or a pair of 
conjugated poles on the second sheet of the complex s plane. Meanwhile 
the simple S matrices do not contain cut contributions from those dispersion 
integrals in Eq. (|2*Tjl. It is not difficult to find solutions of these simple S 
matrices by solving Eq. (J25j) and the solution for each kind of S matrix is 
unique: 

1. A virtual state pole at so with so real. The solution to the scattering 
amplitude is, 



Re R T(s) = 



a V (s R - s ) (s - s L ) 
(sr - s L ) - s ) 



T TY ^ S - Sr) [s L - Sp) 

Im R T S = —— r . (26) 

p(s) {Sr - S L ) [S - So) 



Consequently, the scattering length is 



/ x 2 V /J^ s -s L 

a(s ) = — v - J , (27) 

sr - s L V Sr - Sq 



and the S matrix can be expressed as, 



gfl— SL 



Sis) = ' ^ lizfM^fgzfl . (28) 
l-ip(s)-^- 



s-s L V s R -s 



If there is a bound state at sq, then RcrT in Eq. (|26|) change sign. 
Besides, a bound or virtual state exists only when sl < sq < sr. 
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2. A pair of resonances at z$ (having positive imaginary part) and Zq\ The 
solution for scattering amplitude is, 

Re R T(s) = A( 2o )Ue[^o)]-' S " : 



(s-Zb)(s-2o) 



WT( S ) = , (29) 

and the scattering length is, 

( \ A l \n r c ( M2 (^o) - s fl ) 

a(z ) = A(^ )Re [z p{z )\-f — — - , (30) 

(s R - z ) (s R - z ) 

where 

Im[z Q ] 



A (2b) 



Im[2; ] 2 + Re^op(^o)] 2 
M 2 (z ) = Re[z ] + 



2, . , n..r. 1 , 1,11 M^oP^o)] 



A(z ) = -2Re[^ ] + 



Re[z p(z )} 
lm[zlp{z Q ) 2 } 



lm[z ) 



The S matrix can be expressed as: 



_ M 2 {z )-s + ip{s)sG 
b[ >~ AP(zo)-s-ip(s)sG ' 1 ' 



where 



q = Im [^o] ^33) 
Re[z p(z )) 

Analysis reveals interesting properties of M 2 as a function of Re[^o] for fixed 
Im[zo], as shown in fig. El The inclined straight line corresponds to M 2 = Re[zo] 
which is the asymptotic line of M 2 (zq) when Re[zo] —> oo. The vertical line 
corresponds to Re[zo] = {s R + sl) /2 and it is another critical line on which the 
phase shift 5(s) from a pair of resonances is tc/2 when s — > oo. On the right 
hand side of the line, with the increase of s, the phase shift of the resonances 
can get larger than tc/2 whereas on left hand side of the line, the phase shift of 
resonances can never reach ir/2. Fig. 0] also gives some examples of narrow and 
broad resonances and their contributions to the phase shift. Each resonance 
gives phase shift a positive contribution and the phase shift increases when s 
increases. A narrow resonance will give phase shift a sharp rise and behaves 
much like an ordinary Breit-Wigner resonance, but a broad resonance can only 
give phase shift a slow rise. 
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Re[z ]=(s R +s L )/2 


/' 


-o- 




1 8 2.5- 


i 

i 


E 

0- 




\ 2.0- 
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M 2 (z )=Re[z ] 


i 
i 


-5- 
- J - 


-1.0 -0.5 0.0 


m[z ]=0.5 

0.5- 
0.0- 

0.5 1.0 1.5 2.0 


^ - ' / 



z =1.5+0.1i 
z =1.5+0.5i 

- ' 'z =1.5+1.5i 



Re[z ] s(GeV) 



Figure 4: The left figure shows M 2 ( function of Re[z ], for two choices 

of Im[zo]- The right figure gives some examples of resonances and their con- 
tribution to phase shift. 



4.2 A Pedagogical Analysis to a Toy Model Description 
of Resonances 

A very simple but frequently used parametrization form of S matrix to fit the 
position of a physical resonance is the following, 

where p(s) is the kinematic factor. For equal mass like the nir scatterings, 5* 
matrix as described by Eq. (fH4*j) usually contains a virtual state and a pair of 
resonance poles on the second sheet. But Eq. (J3"4~j) usually contains two pairs 
of resonance poles for nK scattering due to the mass difference between 7r and 
K. Fig. El shows the traces of two pairs of resonance poles with the increase of 
g for different M 2 . We find that M 2 = (s R + s L ) and M 2 = (s R + s L )/2 are 
two critical points for the two pairs of resonance poles. When g increases from 
zero, two pairs of resonances will appear, one from the origin (due to the 1/s 
singularity hidden in the kinematic factor) and another from (M 2 , 0). When 
g increases to a certain magnitude, one pair of resonances will reach the real 
axis and change to two virtual states, and one of them will run upwards s R 
whereas another will run downwards sl- If M 2 > (s R + sl), the left pair of 
resonance poles generated from the origin will change into two virtual states 
and the right pair of resonance poles will become wider and wider but never 
reach the critical line M 2 = (sr + sz,)/2, and if M 2 < (sr + sl) the way of 
motion of the two pairs of resonance poles chang GS, clS shown in Fig. |3| 

Eq. (J34j) is a commonly used parametrization form of 5* matrix to fit res- 
onances. But from the above discussion, we know that Eq. ()34|) is not a 
parametrization of S matrix for one pair of resonance poles, but usually for two 
pairs. Since one pair of poles is below threshold sr or change into two virtual 
states on the real axis (also below the threshold), the existence of such poles 
will violate the validity of chiral expansions at low energies and is therefore 
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-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Re[z ] Re[z Q ] 

Figure 5: The traces of two pairs of resonances from Eq. ()34j) with the increase 
of g for different M 2 . We give two typical figures: the left figure for M 2 > 
(sr + sl), and the right figure for M 2 < (sr + Sl). 



spurious. Actually the low lying poles in here has the same origin as the single 
virtual state pole in the case of equal mass scattering. [IB] They are both gen- 
erated from the kinematical singularity of p(s). Whether the existence of such 
spurious poles strongly affects the determination of another pair of resonance 
poles depends on their contribution to phase shift or scattering length. For ex- 
ample if we set the physical resonance pole position at y/s = 0.8±0.3iGeV then 
another pair of poles will locate at \fs = 0.431 ± 0.125zGeV. The scattering 
length contributed by the two are: a = 2.18GeV~ 1 for k and a = 4.30GeV _1 
for the spurious pole! 

At this moment the major reason to exclude Eq. (J34"|) is rather academic. 
However, in sec. 15.31 we will see that the experimental data also exclude very 
likely the use of Eq. flBH). 

4.3 The Factorized S Matrix and the Separable Singu- 
larities 

For a general S matrix, there can be many poles on the complex s plane. How- 
ever we can always express the complicated physical S matrix into a product 
of many simple S matrices: 

gphy = gcut .JJSf , (35) 

i 

suppose we can find all the poles of S phy . The pole contributions, Sf, to S phy 
are parameterized using the forms given in Sec. 14.11 Now the only uncertainty 
remains is the cut contribution S cut , which, by construction, contains no poles 
but inherits all the cut structure of the original S phy . Hence S cut can be 
parameterized as the following: 

S cut = exp[2tp(s)f(s)] , (36) 
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where f(s) satisfies the following dispersion relation: 

/(.) = fM + ^f -, i>s \) m \ dz 

2,7X1 Jl{z — s)[z — So) 

+ ^-f . lm f Z) r dz, (37) 

7T JR{Z — S){Z — S ) 

where L denotes the left hand cut on the real axis and also the circular cut, 
and R denotes the right hand cut starting from the second physical threshold 
to oo. Apparently the above parametrization automatically guarantees single 
channel unitarity. Furthermore we get, 

disc/ = disc{-i— log \S ph y(s)} } (38) 
2ip{s) L J 

on both L and R. The Eq. f)38j) is derived, from an important property of Sf, 
that is the latter does not contribute to any discontinuity of f. The reason 
follows: firstly we have 



f( s ) = r^log 

J V ; 2ip(s) 



S phy (s)/l[Sf 



J log^(s)- ^£log Sf(s) . (39) 



2ip(s) 2ip 

The first equality in the above equation ensures that / contains no more singu- 
larity than those cuts S phy contains, since poles and zeros of the two S matrices 
on the r.h.s. of the first equality exactly cancel, by definition. From the second 
equality, one easily understands, by comparing with the expressions of Sf in 
sec. 14.11 that the Sf contribution to disc/ vanishes everywhere on all cuts. 
Therefore we have Eq. (J38j) . Especially, when s is on the real axis, Eq. (J38)) 
can be rewritten as, 

l mLtR f(s) = --^-log\S ph y(s)\ , (40) 
2p{s) 

which is a consequence of real analyticity. On the right hand cut R = [(Mk + 
M v ) 2 , oo), Eq. (|40|) can be further rewritten as an analytic expression, 

lm « m = -2W) losvis) = ~^pjs) logi ^i ] ' (41) 

where the subscripts 1 and 2 mean channel 7iK and channel rjK, respectively 
(of course Su = S phy , and in Eq. (J41j) S means the 2x25* matrix of couple 
channel scatterings). In principle, the Eq. (|35|) works not only in the single 
channel region but also works in the inelastic region. However, it should be 
emphasized that all poles in our formulae are on the second sheet. On the other 
side, the phase shift above the second physical threshold is mainly influenced 
by the third sheet pole rather than the second sheet pole (this is true at 
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least for narrow resonances). Our present scheme suffers from the lacking of 
both theoretical 4 and experimental knowledge on the inelasticity parameter, 
77. It is not difficult to imagine the worst situation one may encounter when 
we have only insufficient information on inelasticity: suppose in the inelastic 
region under concern there is no second sheet pole but only a third sheet 
pole. The latter however only shows its effects through Eq. ([41 p. If we neglect 
the inelasticity effects due to our ignorance on rj, we would have to fit the 
phase shift data 5 2 using a second sheet pole. But this is certainly wrong 
by the assumption! Fortunately, the situation just described is not possible 
to happen in ttK scatterings. Because the Kr\ cut is very weak up to the 
Krf threshold [12], the third sheet pole and the second sheet pole should co- 
exist and the mass and the width of the two poles should be similar either. 5 
Therefore, for the if* (1430) pole we are going to introduce in the later fit we 
bear in mind that there are some ambiguities associate with it as discussed 
above but the problem should not be serious. In the later fit we will also 
neglect the Krj' threshold effects which gives some further uncertainties to our 
final results. But the uncertainties should not be large either because we only 
fit the data below the Krj' threshold where the influence from high energy 
cuts should be small and smooth. The uncertainties related to the right hand 
cut and the if* (1430), furthermore, should not waver our main conclusions on 
the k resonance, which is of our major concern in this paper, since the energy 
region where all these problems occur, is rather far from the low energy region 
where the k pole locates. 

Going back to Eq. (|3*Bj) again, it factorizes different singularities of the scat- 
tering S matrix. Also, different contributions to the phase shift are additive: 



where 



i 

p(s)sG[z ] 



6 P = Atan 



M2[z ]-s 

for a resonance located at z (and Zq) and 



(42) 



(43) 



5 P = Atan 



p(s)s 



ls - 


SL 


sr - 


s _ 



(44) 



for a virtual state located at so (sl < so < sr). The bound state contribution 
can be obtained by simply change the sign of r.h.s. of Eq. (j44)l . For the 
background contribution we have 



>BG 



P{s)f{s) 



(45) 



4 In this paper we do not attempt to study a true couple channel problem by giving an 
appropriate parametrization form of rj. The problem was investigated in Ref. |24) but was 
not very successful yet. 

5 This may be best illustrated by a couple channel Flatte model. 
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Notice that the separation of pole contribution and background contribution is 
only a matter of convention. However, our definition of poles and background 
contribution has the advantage - as manifested in Eq. (|40|) - that it greatly 
simplifies the calculations on various cuts. That will be elaborated in the next 
section. 



The approximation scheme in evaluating the cut contributions is to ap- 
proximate S phy in Eq. flUJ) by S xPT on L, 



Since the region where the above discontinuities are being estimated are far 
from the resonance region, we expect the chiral expansion (and hence our ap- 
proximation scheme) works reasonably well here, at moderately low energies. 
Also the logarithmic form of the cut contribution automatically regulates, and 
hence reduces the effect of, the bad high energy behavior of the chiral am- 
plitudes. Actually the cut integrals only need one subtraction (Various cut 
integrals are all subtracted at the nK threshold throughout this paper). From 
these facts, we argue that our approximation scheme for evaluating cuts are 
reasonable, at least qualitatively. The quality of the approximation can be 
tested experimentally, as is done in the next section. 

5 Fits to ttK Scattering Processes 

Our goal is to search for resonances and to study their properties, in ttK 
scatterings, by analyzing the phase shift data. For this purpose, it is very 
important to first estimate various background contributions, which are orig- 
inated from various cut integrals and the subtraction constant appeared in 
Eq. and (|4*2j) . The correct understanding to the background contributions 
is of course greatly helpful in the attempt to a precise determination to the 
pole parameters. It is even vital in answering the question whether there exists 
a broad resonance, like the k resonance under debate. Since the contribution 
from the background and the contribution from the broad resonance can be 
rather similar, the lacking of the knowledge on the former can even lead to 
completely misleading prediction to the latter. It is therefore necessary to 
study carefully the contribution of each term in Eq. ([37)1 . One of the main dif- 
ferences between the contribution from the background and the contribution 
from a broad resonance is their contribution to the scattering length. As we 

1 /2 

will show later that to evaluate the scattering length parameter a generated 
from fits is very helpful in clarifying several issues, including the important 
question whether k exists or not. 

5.1 The Estimation to the Background Contributions 

Our method of estimating the 'left hand cut' contribution is to substitute S phy 
in Eq. (|4H|) by S xPT . The background contribution from the left hand cut is 




(46) 
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then evaluated by calculating the left hand cut integral in Eq. (|S7|) . The once 
subtracted integral is convergent and the integration is formally performed 
from — oo to sl on the real s axis. However, since there is no reliable method 
to estimate correctly the contribution from large negative s region, it is more 
appropriate to truncate the integration along the negative real s axis at A|. 
This introduces an additional cutoff parameter and we will test the dependence 
of the fit results on this parameter. 

The discussions made above clearly shows how approximations enter into 
our scheme. It is necessary to first justify our approximation scheme. Fortu- 
nately it is possible to test whether it is reasonable to use %PT to calculate 
the background contributions. Because in the 1=3/2 channel there are only 
background contributions. So we can compare directly the background contri- 
butions from our approximation with experiments. We have already demon- 
strated in sec. El that the Pade approximation gives poor description in this 
channel to the background contribution. Here we compare Im^,/(s) both using 
S xPT and the [1,1] Pade S matrix, as shown in figs. |U]and[7| 6 We find once 
again that the Pade approximation result is poorer in reproducing the phase 
shift data comparing with the xPT result after adjusting the cutoff param- 
eter within reasonable ranges. Notice that throughout this paper we do not 
attempt to make the global fit by varying those L\ parameters since the L l T 
parameters are in principle threshold parameters. However, several different 
results on the low energy constants found in the literature are tested and it is 
found that the influences from the different choices to the final fit results are 
very small. From fig. [7| disagreement between the theoretical calculation based 
on xPT and the experimental data is found. Better agreement can be achieved 
only by allowing a^j 2 to be much larger than its %PT value. Nevertheless it is 
verified that the problem in the 1=3/2 channel has rather small influence to 
the k pole problem. Further discussions on this point will be given later. 

In the 1=1/2 channel a direct check on the background contribution is 
impossible since this is the channel where we will test the existence of the k 
resonance. However from the experience in the 1=3/2 channel we expect that 
the %PT results on Irri£,/(s) works also well, at least qualitatively. It is worth 
pointing out that the results from the [1,1] Pade approximant are rather similar 
to the xPT results in the 1=1/2 channel, except at large negative s region. 
Therefore we will use the xPT prediction on lmif(s). In the 1=1/ 2 channel, 
there exists also the right hand cut from s = [Mk + M v ) 2 but it is very weak 
till another threshold s = (Mk + M^) 2 is reached. For the right hand cut, %PT 

6 Various theoretical estimations on S^ 2 in fig. [7] are calculated using the following for- 
mula: 



where disci J f(z) are estimated using the tree level, 1 loop x?T results and the Pade result. 
The Adler zero sa for the Pade result is taken as the tree level %PT result on sa- Notice 
that in this way the spurious poles' contributions to the phase shift are absent. 
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Figure 6: Various results on Imx/(s) or / from %PT and Pade approximation 
in 1=3/2 s wave ttK scattering. Left: outside the circular cut; right: on 
L = (-00, -(M| - Ml)}. The L\ constants are taken from Eq. 




! =1.5) 



X PT ,J0 ° P (A 2 =7) 



Figure 7: Theoretical results on 5 versus the data [Sj. The dashed line 
corresponds to the tree level result. Other lines are distinguished from 
the labels in the figure. See footnote El for more explanations. 
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Figure 8: Various results of / and Im^/(s) from ;\PT and Pade approx- 
imation in the 1=1/2 channel. Left: outside the circular cut; right: on 
L = (-00, -(Mi -M$)\. 



result is totally misleading because it violates single channel unitarity in the 
physical region. On the other side, the Pade result maintains single channel 
unitarity but may give a too large contribution to the inelasticity parameter, 
starting from [Mk + M v ) 2 . Since experiments suggest that the inelasticity 
is very small till (Mk + M v >) 2 and the latter is already very far from the 
low energy region we are interested in, we during the fit simply set the right 
hand cut contribution being vanish. As is already discussed in sec. 14.31 this 
approximation may have some effects to the determination of the _fT(1430) 
resonance but should have negligible influence to the k resonance. 

5.2 The Numerical Analysis and Data Handling 

The two experiments on Kir scatterings we will make use of are from the 
LASS Collaboration jl] and Estabrooks et al. [5]. The latter gives the phase 
shift data in the 1=1/2 channel for = 0.73 — 1.3GeV and in the 1=3/2 
channel for yfs = 0.73 — 1.72GeV. It is found that the scattering is purely 
elastic up to 1.3GeV, i.e., r^J 2 = r]^ 2 = 1 in this region. In the 1=3/2 channel 
the inelasticity can be neglected in the whole energy range. The experiment 
by the LASS Collaboration only measures the K~n + channel and hence only 
affords the following combination of data: 

A = aoe* = Tl' 2 + \t^ 2 = - 1) + \^T^ - 1) , (47) 

for y/s = 0.825 — 2.52GeV. Notice that in the above equation the definition of 
T matrix is different from previously used in this paper. 

In this work our discussion will be confined to the single channel approxi- 
mation, though it is understood that our formalism in principle works also in 
the inelastic region. The validity domain of single channel approximation is 
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largely expanded in 7r.fr scatterings since the second threshold, the if 77 channel 
opens very weakly in agreement with SU(3) / expectations, and so any inelas- 
ticity can be neglected until one reaches the Krj threshold ^2]- Therefore we 
in the fit assume elasticity up to the Krj' threshold both in 1=1/2 and 3/2 
channels. As illustrated in the Introduction the main purpose of this paper 
is to study whether there exists the k resonance and, if exists, its properties. 
The approximation of neglecting the inelasticity effects will mainly affect the 
well established if* (1430) resonance which is not our main interest here. 

Our strategy of making the fit follows: assuming two resonances, one for 
if* (1430), another one for the k resonance under investigation and the two 
complex poles contribute four parameters. There are another two parameters 
coming from the two scattering lengths, i.e., aj 2 and a^ 2 (or equivalently, 
two subtraction constants). Therefore there are totally 6 parameters in the fit. 
The first fit only make use of the LASS data up to yfs = 1.43GeV, which is 
about 20MeV below the Krj' threshold and consist of 60 data points. We call 
it Method I hereafter. In the second fit we also include data from Estabrooks 
et al., up to y/s = 1.3GeV in the 1=1/2 channel and up to \fs = 1.43 in the 
1=3/2 channel. This will add another 41 data points, among them 24 come 

1/2 3/2 

from the 5 phase shift data and 17 from the o phase shift data. We call 
this fit scheme Method II hereafter. Treatment to various cuts are already 
illustrated in sec. 15.11 

5.3 The fit to the LASS data (Method I) 

This part of discussion is divided into two subsections: one is the fit without 
considering further constraints from chiral perturbation theory except when 
estimating the left hand cut contributions. The second subsection respects the 
xPT results, especially its predictions on the scattering length parameters. 
Since what is involved here is the SU(3) version of %PT and since there exists 
possible conflict between theory and experiments, we think it is worthwhile to 
be cautious to make the separate discussions. 

5.3.1 The fit without constraints from %PT 

Using MINUIT, we perform the six parameters fit (aj , Q^ 2 and four pole 
parameters for k and if* (1430)) to the LASS data up to 1430MeV. Except 
for these 6 fit parameters we have, as already stressed, one additional but 
somewhat unpleasant parameter: the cutoff parameters responsible for the 
truncation of left hand integrals in both 1=1/2 and 3/2 channels, denoted as 
A|. 7 The fit results will of course depend on the cutoff parameter, nevertheless 
the dependence is much weaker than that of Ref. jTH| which is of course pro- 
gressive. In table El we list various results obtained by varying A| from which 

7 We have also carefully checked the situation when K\ in the 1=1/2 channel and in the 
1=3/2 channel are different. The final results are very similar to those presented in this 
paper. 
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3 


43.55 


0.271 ±0.129 


-0.102 ±0.006 


629 ±81 


697 ± 354 


3.5 


45.96 


0.268 ±0.077 


-0.096 ±0.006 


636 ± 40 


693 ± 236 


5 


53.18 


0.262 ±0.068 


-0.082 ±0.006 


652 ± 30 


685 ± 191 


10 


71.40 


0.254 ±0.062 


-0.060 ±0.006 


676 ± 24 


673 ± 148 


oo 


120 


0.242 ±0.067 


-0.026 ±0.007 


716 ± 16 


653 ± 136 



Table 2: Various fit results to the LASS data up to 1.43GeV obtained by 
varying the cutoff parameters of the left hand integrals. There are totally 60 
data points and 6 parameters. Various results on i^*(1430) are similar and 
are not listed here. All values of mass parameters in the table are in units of 
MeV. 



we draw the following conclusions: 

1. The overall x 2 is rather stable against most changes of the cutoff pa- 
rameters, and the fit prefers A 2 L ~ 1.5GeV 2 . This is supported by an 
independent analysis to the phase shift data provided by Ref. j^J. It 
should be emphasized that these conclusions are obtained only when 
a)J 2 and especially a^ 2 are treated as completely free in the fit. If a^ 2 
is confined to its xPT value, the best A| value will be enhanced. 

2. The fit result on a^ 2 is rather sensitive to the cutoff parameter. Ac- 
tually these sensitivity are only related to the cutoff parameter in the 
1=3/2 channel rather than that in the 1=1/2 channel (see also footnoted 
and fig. El- Remember that we are fitting the combined data now, there- 
fore our method can somehow rather clearly distinguish different channel 
contributions to the LASS data. 

3. Except for the problem mentioned above, most other results are rather 
stable against the variation of A| in a reasonable range. This is even 
true when we confine a^ 2 and d^ 2 in the fit. See later text for more 
discussions. 

To save the already lengthy discussions we in the following will only exhibit 
results for the fixed value of the cutoff parameters: A| = 1.5GeV 2 , unless 
otherwise stated. However we also carefully analyzed the uncertainties of all 
our major outputs induced by the uncertainties of the cutoff parameters and 
found out that the uncertainties are not magnificent. 

As is shown in table El the global fit prefers larger magnitudes of the 
scattering lengths than those predicted by 0(p 4 ) %PT j!3j . which are: a l J 2 = 
0.18 ± 0.02, a ' = -0.05 ± 0.02. The present fit gives however rather large 
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Figure 9: The dependence of the mass, width and their error bars of the k 
resonance on a]j 2 (Method I). The fits contain only 5 parameters. Taking 



A! 



1.5GeV 2 . 



error bars on a l J 2 . To further test the dependence of the k pole on we 

1/2 

also perform the fits by varying a while keeping it fixed during each fit. The 
results are shown in fig. El In fig. El we have not shown the situation when 
a l J 2 > 0.4 where the results become more and more unstable. From fig. El 

1/2 

we may have the impression that when a is roughly less than 0.35 the k 

1/2 

resonance seems to exist. If a is greater than the value the existence of the 
resonance becomes doubtful. This observation will be further examined in the 
following discussions. 

The full results of the six parameters fit to the LASS data are the following: 



xLj. = 38.35/(60- 6) ; 
M K = 594 ± 79MeV , T K 



1/2 



0.284 ±0.089 



3/2 



M K * = 1456 ± 8MeV , Y K 



724 ± 332MeV ; 
= -0.129 ±0.006 ; 
= 217±31MeV r . 



(4? 



The corresponding fit results are plotted in fig. Since their is no separate 
data of 1=1/2 and 1=3/2 channels provided by the LASS Collaboration, we 
also plot curves of and S^ 2 from our fit results versus the phase shift data 
from Estabrooks et al. in fig. ^2 We stress again that in the present fit the 
LASS data are combined data from both the 1=1/2 and 1=3/2 channels. From 
fig- EH we find that the present method can rather clearly distinguish different 
contributions from different channels. 

In order to answer the question whether k exists or not, we further freeze 
the k degrees of freedom in the fit and we get 



XL./. = 63.67/(60-4) 
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Figure 10: Fit results (Eq. (|1HJ)) to (fio (left) and a (right) using Method I. 
The data are fitted only up to 1430MeV but are shown up to 1500MeV. 
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M K = , T K = ; 

a Q /2 = 0.446 ± 0.006 , a^ 2 = -0.130 ± 0.006 ; 

M K * = 1432 ± 3MeV , T K * = 314 ± 20MeV . (49) 

Comparing with the results in Eq. (|4*Hj) the Xd.o.f. given by Eq. (j4*9~j) is increased 
by a factor of 1.7. If this is not enough to support the existence of k, the value 
of a l Q 2 given by Eq. (J4T?|) is too large comparing with the xPT value. While the 
fit value of a?J 2 is about 4 a away from %PT result, predicted by Eq. (f4~9"|) 

is 14cr away! It is actually easy to understand why the fit program is forced 
to chose such a large a^ 2 . Looking at fig. when the contribution of the 
k resonance is withdrawn the subtraction constant has to be much increased 
to fit data. Since the pfo contribution increase faster at threshold than the k 
contribution and slower at higher energies than the k contribution, this leads 
to a larger a]j 2 . 

It is also worthwhile to investigate the fit using Eq. ()34j) to parameterize 
k. The result follows: 

XL./. = 38.22/(60- 6) ; 

M R = 726 ± 8MeV , T K = 698 ± 19MeV ; 

a /2 = 0.403 ± 0.025 , aj 2 = -0.129 ± 0.006 ; 

M K * = 1453 ± 3MeV , IV = 220 ± 18MeV , (50) 

and the spurious resonance pole locates at (the central value): M = 477MeV, 
T = 94MeV. Comparing with Eq. (}4~%)) . the total \ 2 given by Eq. (joT?|) is 

1/2 

very similar, but it gives a too large a comparing with the %PT prediction. 
Therefore Eq. (jH4^) is very likely to be excluded by the LASS data. 

5.3.2 The fit with constraints from %PT 

It is sometimes found in the literature the discussions on the constraint of the 
Adler zero on the scattering amplitude. [23 In our scheme, it is possible to 
embed this constraint into our parametrization form, similar to what is done 
in Ref. ^7j. However, the Adler zero automatically emerges in our approach 
if the subtraction constant fo is limited within certain range. This is because 
all the resonance(and virtual state) S matrices are real and less than 1 when 
Sl < s < sr. On the contrary the cut integrals contribute a factor larger than 
1, therefore a T matrix zero emerges in the right place when f is confined to 

a certain range, which in turn put some constraints on the magnitude of the 

l li 

scattering length parameter itself. For example, within the range < a > < 
0.26 (see fig. EJ) there exists a T matrix zero in the region sl < s < sr, and 

1/2 

when a ~ 0.20 the zero locates in the place close to the one loop %PT 

1/2 

prediction. We make further fit by confining a in the region 0.18 ± 0.02 and 
the results follow: 

XL./. = 38.96/(60- 6) ; 
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M K = 646 ± 7MeV , T K = 540 ± A2MeV ; 
a Q /2 = 0.2 , al /2 = -0.128 ± 0.006 ; 

M K * = 1450 ± 5MeV , T K , = 232 ± 25MeV . (51) 

The Adler zero position is now at sa — 0.245GeV to be compared with the 
1-loop xPT value sa — 0.233GeV. The necessity for the existence of the 

1 /2 

k resonance may be best illustrated by the fit when constraining a to lie 

1/2 

within the range as predicted by value, a = 0.18 ±0.02 and meanwhile 
freeze out n. Under this situation a[j will reach its upper value at 0.2 and the 
Xtot ~ 750! Comparing with the value of Xtot m Eqs. (}5T|) and (J4*8|) . it clear 

1 /2 

demonstrates the necessity to include the k resonance if a is close to (or not 
much larger than) its xPT value. 

The next question is to ask what would happen if we further confine both 
a\j 2 and a^ 2 to their x?T value? The answer is, 

XL./. = 127.6/(60- 6) ; 

M K = 655 ± 9MeV , T K = 549 ± 42MeV ; 

a /2 = 0.2 , a /2 = -0.07 ; 

M K * = 1465 ± 5MeV , T K * = 258 ± 33MeV . (52) 

Though the Xtot ^ s now m uch enhanced as comparing with the result given by 
Eq. (joTjl . we are very much consoled by the nice agreement on the k pole be- 
tween the two methods. This again suggests that though we fit the LASS data 
which is a combined effect of both the 1=1/2 channel and the 1=3/2 channel, 
the present approach can somehow clearly distinguish different contributions 
from different channels. Actually if we make a plot like fig. ^2 we find that the 
fit curve in the 1=1/2 channel still agrees well with the data whereas in the 
1=3/2 channel there exists rather large deviations. Technically, the large Xtot 
appeared in Eq. ()52|) can be reduced roughly by half when increasing K 2 L (or 
more precisely, increasing the cutoff parameter in the 1=3/2 channel). This 
ambiguity with respect to the cutoff parameter will contribute, though not 
large, some uncertainties to the pole positions, which can be considered as the 
'systematic' error in our approach. 

5.4 The combined data fit (Method II) 

The data from Estabrooks et al. contain some problems. This can be clearly 
seen from fig. [TT] The dip structure of the phase shift data around 920MeV 
apparently leads to j- < 0. According to our analysis, it cannot be explained 
by the smooth background contributions. If the dip structure were true it 
unavoidably leads to very complicated pole structures and causality violation 
conclusion. The dip structure actually gives a large contribution to the total 
X 2 in the combined fit to both the LASS data and the data from Estabrooks et 
al.. Nevertheless we still make the combined fit using the method as described 
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Figure 12: The combined fit results (Eq. ([53))) on 5 t 
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previously. The result for six free parameters fit follows: 



X 



(I 



1/2 







M, 



K'' 



266.9/(101 - 6) ; 
= 898 ± 315MeV 
= 0.384 ±0.017 , 

= 1468 ± iMeV 



3/2 



K* 



1902 ± 282MeV ; 
-0.136 ±0.002 ; 
= 146 ± 23MeV . 



(53) 



We see from the above equation that the Xdof * s m uch larger than the fit 
using only the LASS data and the scattering length parameter is much larger 
than the %PT value and the result should not be trustworthy. 8 We also plot 
the fit results versus the 5q 2 data in fig. El It is interesting to notice that the 

1/2 

results of M K and T K and their error bars obtained by varying a , as shown 
in fig. have a very similar behavior to the results obtained using Method 
I as shown in fig. EJ for small value of a)J 2 . Therefore, though not successful, 

1/2 

the combined fit still confirms the conclusion that the k resonance exists if a 
does not deviate too much from its ;^PT value. 



6 Summary and Conclusions 

We have already made a rather long and detailed discussion on the new uni- 
tarization approach and the results on the k pole based on the approach. In 
here we summarize our main physical results: first of all, there exists the k 
resonance, if the scattering length parameter a]J 2 does not deviate much from 
its xPT prediction. Secondly, according to our fit to the LASS data, the mass 
of the k resonance is definitely smaller than most previous results found in the 
literature. It is found that the width parameter is more flexible than the mass 
parameter in the fit results, which can also be seen in fig. El Our fit results are 

furthermore, for the large fit result a^ 2 there is no Adler zero. 
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Figure 13: The dependence of the mass, width and their error bars of the k 
' (Method II). Taking Af = 1.5GeV 2 . 
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resonance on a 



given in Eq. (|48j) . If we further fix a to be its current xPT value within la 
error bar, the pole position of k is approximately given by Eq. (J5T|) . However, 
unlike tttt scatterings, the scattering lengths under concern are from SU(3) 
%PT results and one expects that high order corrections are more significant 
here. Therefore one should also be cautious when referring to Eq. (|5ip or (|52|). 
Furthermore, there are some mismatch or disagreement between our method 
and numerical results and the experimental data which predominately come 
from the 1=3/2 channel. 9 However, we have made efforts to separate and to 
minimize the ambiguity when fixing the location of the k pole. In fact, we 
are convinced by comparing Eqs. (JoTj) and (p)2*|) that the problem in the 1=3/2 
channel has only minor influence to the K pole location. 

Our numerical results, especially Eq. (p)2~j) and the way they are derived 
are in principle consistent with that of Ref. ^2] m which k pole with similar 
location can be found after adding a few 'data points' generated by %PT to 
the original LASS data. These 'data points' should take the similar role as 
constraining the scattering lengths here. The present approach also share some 
similarities, at qualitative level, with that of Ref. [7] , though in details the two 
approaches differ in the treatment on both the pole contributions and the cut 
contributions. 

While this paper is being completed, we received a paper by Biittiker, 
Descotes and Moussallam |27j . The central value of a l J 2 given in Ref. [2Z! 
obtained by solving Roy-Steiner equation increases considerably than that of 
Ref. [13J, and a^ 2 changes very little. Apparent disagreement on 5^ 2 between 
the result of Ref. and the experimental data are also indicated. If we 

9 The discrepancy may be reduced in future when a two loop xPT calculation becomes 
available. This expectation is supported by a similar study on tttt scatterings 26 . 
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7 


69.8 


0.218 ±0.026 


-0.053(at limit) 


682 ± 25 


598 ± 153 


9 


71.1 


0.234 ±0.032 


-0.053(at limit) 


682 ± 24 


629 ± 134 


11 


74.6 


0.243 ±0.032 


-0.053 ±0.009 


683 ± 16 


651 ±85 


20 


90.1 


0.245 ±0.021 


-0.044 ±0.006 


693 ± 14 


658 ± 40 


oo 


137.4 


0.245(at limit) 


-0.038(at limit) 


710 ±14 


661 ±29 



Table 3: Various fit results to the LASS data up to 1.43GeV obtained by 
varying the cutoff parameters of the left hand integrals and with the constraints 
on a^ 2 and a^ 2 from Ref. [2Tj . All values of M K and T K in the table are in 
units of MeV. 



1 /2 

take roughly a ~ 0.201 — 0.245 as indicated by the results of Ref. [27] (also 
d^ 2 is readjusted to lie within the range -0.053 - -0.038), we would obtain, 
similar to obtaining Eq. (|52j) . table El Again, when there is a large x 2 occurs in 
table H3 it is mainly contributed by the 1=3/2 channel. The 'systematic' errors 
induced by varying cutoff parameters are estimated and the final results on the 
k pole position obtained by considering the constraints provided by Ref. [23 
is estimated from table |HJ 

M K = 688 ± 25 ± 22MeV , T K = 613 ± 153 ± 48MeV . (54) 

Here the first error bar corresponds to the largest error bar generated by MI- 
NUIT when varying the cutoff parameters. The second error bar is our estimate 
obtained by examining the variation of the corresponding central value when 
changing the cutoff parameters, which is also a rather conservative estimate. 
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